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I. INTRODUCTION 

The fractional quantum Hall (FQH) effect^ is the 
archetypal system for the emergence of topological or- 
deJ^ in condensed matter physics. Due to the presence 
of strong correlations, the theoretical understanding of 
its microscopic properties heavily relies on finite-size nu- 
merics.l^^ Such calculations are limited to rather small 
system sizes, as the many-body Hilbert space grows ex¬ 
ponentially with the system area. An ingenious way to 
circumvent this problem comes from the recen t develop- 
meniPtSl of exact matrix product stateJ^^^ (MBS) for 
a large class of FQH model wave functions derived from 
conformal held theory (CFT) correlator Exploit¬ 
ing the area law of quantum entanglement,!^ the MBS 
factorization enables efficient calculation of physical ob¬ 
servables and gives access to much larger system sizes 
than previously attainable. 

In a recent paper,^^ we applied this technique to 
study quasihole excitations in the Moore-Read,the 
Gaffnian,!^ and the Z 3 Read-RezayP^ states. These ex¬ 
otic quasiparticles were conjectured to be non-Abelian 
anyons, and they constitute the most striking ma nifest a- 
tion of the topological order in the FQH liquids A 
hallmark of their non-Abelian character lies in the topo¬ 
logical degeneracy of multiquasihole states. For each set 
of hxed quasihole positions, there exists a quasidegener¬ 
ate subspace of states, protected against local perturba¬ 
tions, and braiding the quasiholes induces unitary trans¬ 
formations over this degenerate subspace. These trans¬ 
formations depend only on the topology of the braids, 
rather than their actual shapes, and those induced by 
distinct braids do not commute. In Ref. na using the 
MBS representation, we explicitly demonstrated for the 
first time the Fibonacci nature of the Z 3 Read-Rezayi 
quasiholes from a microscopic calculation. We estimated 
the quasihole radii and quantified the length scales asso¬ 
ciated with the exponential convergence of the braiding 
statistics, but did not provide details regarding the con¬ 


struction of the quasihole wavefunctions. In this paper, 
we discuss in detail the novel technical aspects of the 
MBS representation for the non-Abelian quasiholes. 

The construction of the exact MBS is based on the 
rewritin g of F QH model wave functions as conformal cot- 
relators.^ ^^ l ^^ l This elegant formalism provides a particu¬ 
larly nice way to resolve the topological degeneracy of 
non-Abelian quasih olein terms of the so-called con¬ 
formal blocks.!^^^^ Each conformal-block wave function 
is indexed not only by the quasihole positions, but also 
by a tree of topological charge labels specifying the dif¬ 
ferent fusion channels of the quasiholes. Enumerating 
all the possible fusion tree labelings compatible with a 
given theory generates a complete basis over the degen¬ 
erate subspace. A special benefit of the conformal-block 
basis is the explicit manifestation of the putative braid¬ 
ing s tatist ics in the analytic structure of the wave func- 
tions.lsMH Specifically, as a function of the complex quasi¬ 
hole coordinates, the conformal blocks display branch-cut 
singularities emanating from the quasihole centers. The 
monodromy matrix associated with crossing the branch 
cuts is conjecturecP^ to coincide with the corresponding 
quasihole braiding matrix, up to an Abelian Aharonov- 
Bohm phase due to the magnetic field. This conjec¬ 
ture rests on the observation that the overlaps between 
conformal blocks resemble the partition function of a 
classical plasmeP with pinned quasihole charges, a pe¬ 
culiar feature of the quasihole-dependent normalization 
of the conformal blocks.!^ This link eliminates the need 
to directly integrate the non-Abelian Berry connection 
to compute the braiding matrix.f^SHlIl Xo demonstrate 
that the braiding statistics is indeed captured by the 
monodromy of conformal blocks, we only have to estab¬ 
lish (through a microscopic calculation) that the plasma 
is in a screening phase. This simplification led to the 
analytic identihcation (with the assumption of plasma 
screening) of the Moore-Read quasiholes as Ising anyons 
in Ref. 1321 and also played a crucial role in our numeri¬ 
cal demonstratiorfi^ of the Z 3 Read-Rezayi quasiholes as 
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Fibonacci anyons. The main purpose of the current pa¬ 
per is to explain how to translate the conformal blocks 
into calculation-friendly MPS form while preserving the 
two highly desirable features, namely, the monodromy 
structure and the plasma normalization. 

The paper is organized as follows. In Sec. |n) we pro¬ 
vide a pedagogical review on the construction of quantum 
Hall MPS in the absence of non-Abelian quasiholes, and 
in particular, we derive the plasma normalization from 
conformal correlators on the cylinder. And as a precur¬ 
sor, we also discuss the MPS representation of Abelian 
quasiholes. In Sec. in we proceed to the non-Abelian 
case. We explain step by step how to construct MPS 
for the conformal-block wave functions with quasiholes 
in the bulk. A large part of the discussion is devoted to 
the derivation of the subtle commutation rules between 
a non-Abelian quasihole insertion and the electron oper¬ 
ators. We provide explicit recipes for the Moore-Read,l^ 
the Gaffnian,!^, and the Z 3 Read- RezayP^ states. Tech¬ 
nical details of the construction are addressed in the ap¬ 
pendices. 

II. MATRIX PRODUCT STATES FROM 
CONFORMAL CORRELATORS 

Before discussing the non-Abelian quasiholes, we first 
review the constructiorP®^^ of the quantum Hall matrix 
product states (MPS) from conformal correlators. In 
preparation for the later discussion of the quasihole wave 
functions, we pay special attention to preserving the nor¬ 
malization of the conformal correlator. 

We consider model wave functions in the lowest Lan¬ 
dau lewl constructed from chiral conformal correla- 
tors.GsUll jn formalism, an electron at position z is 
represented by a primary field insertion V(z), and the 
conformal correlator 

iv{Zi)V{z2)---V{Zn)) ( 1 ) 

can be viewed as a many-body wave function 
^ 2 , ■ • •,-Zn)- Here and hereafter, the single brack¬ 
ets (•) denote a CFT correlation function, in contrast 
to the double brackets ((•)) representing the states (wave 
functions) of the physical electrons. For example, the 
Laughlin wave functiorP at filling v can be describecP^ 
in terms of a massless free boson (j). The electron operator 
V{z) is the normal-ordered exponential 

V(z) . (2) 

Using the propagator {(j){z)(p{z')) = — log(z — z') in the 
plane, we recover from Eq. Q the familiar Laughlin wave 
function ~ more complicated quan¬ 

tum Hall states (see Sec. |HI[), the corresponding CFT has 
a direct-product structure® where in addition to the free 
boson, we also have a separate so-called “neutral” CFT. 

From now on we adopt the cylinder geometrj® with 
finite perimeter Ly. The complex coordinate z = x -\- iy 


has X running along the cylinder axis and y around its 
perimeter. For convenience, we set the magnetic length 
to unity, and we define the inverse cylinder radius 

27r ^ ^ 

7=i-- (3) 

Ijy 

The many-body wave function 'i'{zi, Z 2 , ■ ■ ■, Zn) is given 
by the conformal correlator in Eq. Q evaluated in the 
cylinder geometry, which can be mapped the usual planar 
geometry through the conformal transformation z —>■ . 

Interpreting the x coordinate as the imaginary time, 

the CFT Hamiltonian is given by 7(^0 ~ ^c), with Lq 
being the Virasoro generator for dilations and c being the 
chiral central charge. For the direct-product theory of a 
neutral CFT and a free boson, we have the decomposition 

f f neut I f boson rA\ 

Lq — Lq Lq . (4j 

In this section, we focus on the boson part. The mode 
expansion of the chiral field (p on the cylinder is given by 

(l){z) = $o-ilzao + iy^^- ( 5 ) 

n^O 

Here, the a„ modes of the U(l) current satisfy the Heisen¬ 
berg algebra 

[hn, UTn] — Ti (6) 

while 4>q is the canonical conjugate to the zero mode oq, 
[<^o,ao]=L (7) 

And the dilation operator for the free boson is given by 

^boson ^ ^ ^ 1-2^ (8) 

m>0 

The ao operator measures the U(l) charge in unit of 
times the electron charge, in the sense that 

[V(z)]"^ ao V(z) = ao-b T. (9) 

The zero mode operators 0o and ag are decoupled from 
(commute with) the ladder operators a„^o- As a result, 
the free boson Hilbert space can be split into sectors la¬ 
beled by the U(l) charge do, with 0o coupling different 
sectors. The primary state IQ) = :e*‘5‘^°:|l) has U(l) 
charge Q, and the corresponding Hilbert space sector 
with charge Q is spanned by the descendants of |Q) under 
the boson modes {a„<o}. 

A. Background charge and gauge choice 

The MPS is a tensor factorization of the second- 
quantized amplitudes of the many-body wave function 
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'i>{zi,Z 2 ,...,ZN) in Eq. gpi 

The first step of the con¬ 
struction is to obtain in the occupation-number basis. 
We choose the Landau gauge for the magnetic field, and 
work with the orbitals labeled by the wave number j € Z 
in the y direction: 




{V^Ly)h 




( 10 ) 


These one-body states take the form of a holomorphic 
function in 2 ; times a Gaussian in x. Due to the chirality 
of the electron operators V(z), the conformal correlator 
in Eq. Q does not produce the (non-holomorphic) Gaus¬ 
sian factor, and thus does not yet qualify as a many-body 
wave function in the lowest Landau level (in any gauge). 
Fortunately, the Gaussian factor can be generated nat¬ 
urally by spreading the neutralizing background charge 
for the boson field (j) uniforml'^^ on the cylinder. This 
amounts to inserting another (non-primary) field 



into the conformal correlator, representing the neutraliz¬ 
ing background charge at filling v. Here, the integration 
is performed over the cylinder surface, and the normal 
ordering removes unwanted interactions between back¬ 
ground charges at different locations. However, as dis¬ 
cussed in Ref. [m this extra insertion of the background 
charge has a side effect: in addition to the desirable Gaus¬ 
sian factor, it also introduces a non-holomorphic gauge 
factor. Taken altogether, the cylinder many-body wave 
function in the Landau gauge is given bjEl 

4'(zi,...,z„) =e*^-"'*''(v(zi)---V(2„) Obc). (12) 

This relation is proved in Appendix [X| 


B. Occupation-number basis 


• • ,ju} as 




n 


dy^ 

Lii 


'i’{lji + iyi,...,'yjn + iyn) 


= n / ^ + m) ■ ■ ■ V(7jn + iyn) Obc), 

. Jo ^y \ ' 

(14) 

up to a constant normalization factor. Here and hereafter 
we consider only the fermionic quantum Hall states. No¬ 
tice that at Xi = jji, the gauge transformation in 

Eq. (12) cancels the Fourier factor 

The next step is to rewrite the above expression in 
terms of the occupation numbers rrij = 0,1 of each Lan¬ 
dau orbital j. We work in the operator formalism of 
the CFT, and interpret the x direction along the cylin¬ 
der axis as the imaginary time and the perpendicular y 
direction as the space direction. For simplicity, we first 
consider only the electron operators, and postpone the 
treatment of the background charge operator. We are 
free to pick the index ordering of the occupied orbitals. 
Choosing ji > • > jn ensures the time ordering 

of the electron operators and allows us to convert the 
correlator into an operator expression 


(v( 7 ji +iyi)-- -Vijjn + iyn)^ = 

(out|V( 7 ji -k iyi) ■ ■ ■ Vi'jjn + *2/n)|in). (15) 

We use over-hats on the right-hand side to highlight the 
operator nature of the insertions, and we can set the 
in- and the out-states to the vacuum. Thanks to the 
conformal invariance, the x dependence of the primary 
field insertion V{x + iy) can be isolated, 

V{x + iy) = , (16) 


with 7 L 0 being the CFT Hamiltonian on the cylinder. 
We define the zero-mode of the electron operator as 


Vo 





(U) 


We now try to expand the above wave function in 
the occupation-number basis. Since each Landau orbital 
is a momentum eigenstate, we can extract the second- 
quantized amplitudes through a Fourier transform in the 
y direction. Notice that the one-body wave function ipj 
reduces to a simple plane wave along the orbital center 

X = 7j, 


( 13 ) 

Taking advantage of this,l^ we place the Fourier in¬ 
tegration contours along the orbital centers, and ex¬ 
press the amplitude associated with the occupied orbitals 


Without worrying about the background charge for now, 
we can rewrite Eq. (14) as 


(out| Vo Vo ■ ■ ■ 

Volin), (18) 

up to factors that depend only on the energy of the in- 
and the out-state boundaries. The above expression is an 
imaginary time evolution along the cylinder axis, punc¬ 
tured by the electron zero-mode operators at the center 
of each occupied orbital.l^ To assign the time evolution 
to individual orbitals, we define 

U{s) = 


(19) 
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which advances in imaginary time any CFT state by s 
along the cylinder. Finally, we can write the second- 
quantized amplitude (({ni}|'I')) associated with the occu¬ 
pation numbers {m} = [mo, toi, m 2 ,...] of the Landau 
orbitals as 

(({m}|4')) = (out|---C""=C'”^iC'™«|in), ( 20 ) 

with the orbital C"* operators given by 

C° = C/(7), Cl = 17(7) Vo- (21) 


Upon choosing a basis for the CFT Flilbert space, the 
operator expression in Eq. (20) becomes a matrix product 
state. 


C. Cylinder evolution operator 


We now go back to the issue of the uniform background 
charge. We would like to treat it in the same way as the 
electron operator. To this end, we first split the two- 
dimensional integral in Cbc [Eq. 0 ] into small patches, 

J<ixdy${x+iy) . . ^-i^SxSy (ji{x+iy) . ^22) 

x,y 


where the product over patches (x, y) is time ordered, 
and V is the filling fraction. Evidently this opera¬ 
tion introduces unwanted self-interactions between back¬ 
ground charges at different locations. Fortunately, this 
only adds an overall constant factor that does not de¬ 
pend on the electron position. Notice that each fac¬ 
tor : : is now a primary field, to which 


Eq. (16) applies: 


.Q-i^SxSy<j>(x+iy). _ ^X’jLq . ^-i^SxSy <j>(iy) . ^-X’yLg 


Thanks to the time ordering, we can recombine the 
patches at the same x, and up to an overall constant 
we have, 

■ Jdxdy^{x+iy) . ^ ga;7Lo g-* 00 g-a;7Lo ^23) 

X 


where the product over time slices x is still time ordered, 
and the zero-mode of the boson field [Eq. ([^] is picked 
up by 

$0 = (24) 

J-Ly/2 dxy 

Therefore, up to an inconsequential overall constant, 
the insertion of the uniform background charge operator 
amounts to injecting an exponentiated boson zero mode 
at each time slice. We can combine this with the time 
evolution, and redefine U{s) as the path-ordered expo¬ 
nential 


U{s) = V exp 



^ 7-^0 + i 



(25) 


where Lq is the dilation operator [Eq. ([^] for the direct- 
product CFT. This modification is enough to capture 
the effect of the uniform background charge operator. 
Note that in Eq. (25) the path dependence comes from 
the boson zero mode (po and its canonical conjugate oq 
hidden in Lq [Eq. ^]. This allows us to simplify the 
path-ordered exponential,^ yielding 


Uis) = 


.S^/V . , 

exp I —i - 00 exp 


-S7L0 {Vvao + 


(26) 


where the two exponentials do not commute. This 
new expression for U{s) supersedes the original defini¬ 
tion in Eq. (19), and it enters the second-quantized am¬ 
plitude through the orbital operators (7° = [/(y) and 

Cl = 17(7 )Vo [Eq.@]. 

We emphasize that the above treatment of the cylinder 
evolution operator is not specific to correlators of the 
electron operator. The resulting formula for U{s) applies 
generally to time-ordered correlators of any conformal 
primary fields in the presence of a background charge. 
We will make use of this fact when we derive the MPS 
representation of the quasihole insertion in Section II E| 

Recall that oq measures the U(l) charge in unit of -/F 
times the electron charge. Using Eq. 0 , we find that 


[f7(s)]-iao[/(s)=ao-—• (27) 

7 

Letting s be the orbital spacing 7 , we see that the amount 
of background charge contributed by each Landau orbital 
is equal to —v times the electron charge. This indeed 
neutralizes the total electric charge at filling iz. 


D. Matrix product factorization 


The second-quantized amplitude in Eq. (20) can be 


readily converted into a matrix product state. Between 
each pair of adjacent (7™ operators, we can insert a unit 
resolution into a complete set of states over the conformal 
Hilbert space. 


i = (28) 

oc 

For the free boson, the orthonormal basis states |a) are 
simply normalized descendants under the U(l) current. 
For the neutral CFT, the Virasoro descendants are in 
general not orthogonal, or even linearly independent. An 
orthonormal basis can be obtained through the Gram- 
Schmidt process after eliminating the null modes.^ 

Due to the e~'^ ° factor introduced by the cylinder 
evolution [/(y) in the (7™ operators, the CFT states with 
higher energy (as measured by yLp) are exponentially 
suppressed at finite cylinder perimeter. This allows us to 
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(Jx + 1)7 X Jxl 


FIG. 1. Insertion of a single quasihole at imaginary time y. 
The solid vertical lines in orange mark the center positions of 
the two Landau orbitals sandwiching the quasihole. 


truncat^^ the conformal Hilbert space by keeping only 
the lowest few levels of descendants. The resulting finite¬ 
dimensional vector space is the MPS auxiliary space, over 
which the orbital C"* operators assume a matrix repre¬ 
sentation 


[c^U = 


(29) 


The calculation of these matrix elements is discussed 
in Appendix With the truncated representation of 
the C"" operators, the second-quantized amplitude in 
Eq. (20) becomes a product of matrices dotted into the 
boundary vectors, which can be evaluated numerically.^] 


E. Abelian quasiholes 


In the GET formalism,^^ similar to the electrons, a 
localized quasihole at 77 S C is represented by a primary 
field insertion in the conformal correlator. As a warm-up 
for the non-Abelian case, in the following we first discuss 
the Abelian quasihole, represented by 


2 ( 77 ) =:e 




(30) 


This operator couples only to the free boson, and it gen¬ 
erates the familiar quasihole ~ v) factor when eval¬ 

uated in the plane (as opposed to the cylinder). We now 
discuss how to add 2 ( 77 ) to the MPS construction.!^ 

We consider a single Abelian quasihole at r] = x + i(- 
In Eq. (15), we now have an extra 2(x + *C) operator in¬ 
serted into the chain of electron operators, at the position 
determined by time ordering. Similar to Eq. (16), we can 


extract the x dependence using the dilation operator. 


2(x + *C) = e^^'^“2(7C)e 


(31) 


Following the same steps from Sec. |IIB to Sec. |IIC| for 
each imaginary-time interval of size s between adjacent 
primary field (electron or quasihole) insertions, we can 
capture the dilation and the background charge using 
the cylinder evolution operator U{s) [Eq. ([^]. 

As a result of the quasihole at x +J'Cj the cylinder 
evolution is now further punctured by 2(*C) ^-t the time 
slice X- To locate this time slice relative to the Landau 
orbitals, recall that the orbital j is centered at x = qj 
[Eq. (|Io|]. The 2 ( 7 C) operator is thus inserted between 
the orbitals and + 1 (Fig. [^, with 


Jx = belli ■ 


(32) 


Here, the floor function [tj G Z denotes the largest inte¬ 
ger no greater than t G U. The insertion of 2(*C) breaks 
the cylinder evolution I7(j) associated with the orbital 
C"* operator at into two parts, 

mi) ^ mm)m)mi-m), (33) 

where Sx G (0, 7 ] denotes the displacement of x from the 
center of the orbital -I- 1 (Fig. [^, 

<5x = Ox + 1)7 - X- (34) 


Note that the cylinder evolution operator (7 (s) does not 
commute with the QOC) insertion, and from its definition 
in term of the path-ordered exponential in Eq. (25), we 
have 


uii-sx) = [u{sx)]-^ mm 


(35) 


Then, without modifying the C™ operators, the quasi¬ 
hole 2(x + *0 can be represented in Eq. (20) by the 
insertion of 


U(5x) 2«) [U{5x)m (36) 


between the C"" operators for orbitals jx and jx + ^- 
There is one extra complication that we have glossed 
over. When we expand the CFT-derived wave function 


V(zi)---V(z„)2(x + <) Clbc 


(37) 


into Slater-determinant basis states in Eq. (14), we place 


the electron contours at the center of each occupied or¬ 
bital. To bring the field insertions into time ordering, 
the electrons on orbitals with center position x < x have 
to be moved across the pinned quasihole field at time x- 
For the Abelian quasihole as in Eq. (30), each of these 
commutations incurs a minus sign. 


(• • • V(z)2(x+<) •••) = -(••• 2(x + <)V(z) • • •). (38) 


The above anti-commutativity simply reflects the fact 
that in the planar wave function an Abelian quasihole 
is represented by the odd-power factor ~ ‘T)- Tor- 

mally, one could also derive this minus sign from the 
operator product expansion. 

We need to collect these minus signs together and at¬ 
tach them to the quasihole insertion. To this end, for each 
Slater-determinant basis state we have to count the num¬ 
ber of occupied orbitals with center position x < x- This 
number can be extracted from the conserved U(l) charge 
at the time slice x of the quasihole insertion. [Note that 
for each orbital the cylinder evolution 1 /( 7 ) contained 
in the (7™ operator correctly accounts for the associated 
background charge.] Specifically, the number of occupied 
orbitals with x < x is given by 

\/^ao + Ox + 0^ 


(39) 
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inserted between the C™ operators for orbitals and 
+ 1. Here, the zero-mode oq measures the U(l) charge 
in unit of ^ times the electron charge e, and the second 
term cancels the background charge —ve carried by each 
orbital j S {0,1,..., j^-}. Finally, we can write down the 
full expression for the quasihole operator in the MPS 

um Q{<) [UiSx)]-" (_i)v/^ao+a,+i)._ (40) 

Note that og does not commute with as the 

latter contains 4>q. 

The above construction can be easily generalized to 
the case of multiple Abelian quasiholes, with 

• • • V{Zn) Q(Xl+<l) • • • Q(Xm+<m) C>bcy (41) 

We put the m quasiholes in time ordering. 


branch cuts originating from the center of each quasihole. 
We will discuss this in details in Sec. IIIF Finally, we 
emphasize that our prescription for the Abelian quasihole 
matrix differs from Ref. IHl in the more careful handling 
of the background charge. The new formula here exactly 
preserves the quasihole-dependent normalization of the 
conformal correlator. This feature is desirable as it en¬ 
ables us to leverage the plasma analogy when checking 
the braiding statistics.!!^ 


III. NON-ABELIAN QUASIHOLES 

We now proceed to the construction of MPS for the 
non-Abelian quasiholes. We consider the so-called (fc, r)- 
clustered state^^ at filling fraction u of the lowest Lan¬ 
dau level, with k, r being integers and 


Xi > X2 > ■ • • > Xm- (42) 


Then, each quasihole has the MPS representation given 
by Eq. (40), except for one minor modification. For the 


l-th quasihole, to extract the number of occupied orbitals 
with a: < xz I "we need to subtract not only the background 
charge, but also the U(l) charge introduced by the other 
quasiholes inserted on its right. This modifies the com¬ 
mutation sign in Eq. (40) to 




(43) 


The last two terms in the exponent [(j^^ -|- l)z^ — (m —l)v] 
only introduce an overall constant phase factor (inde¬ 
pendent from the electron occupations), but they are 
necessary to eliminate the branch-cut ambiguity in the 
fractional power (—1)^“°. This branch cut ambiguity 
comes from the fact that the operator ^JvaQ now counts 
fractional quasihole charges. 

This finishes our review of the MPS construction in 
the absence of non-Abelian quasiholes. We conclude this 
section with a few remarks. Eirst, the quasihole MPS 
has a rather subtle parametric dependence on the quasi¬ 
hole position x + *C- This dependence is not holomorphic 
due to the presence of the nonchiral background charge 
operator Ohc, and the x-dependence differs considerably 
from (. The imaginary part C enters the MPS directly 
(and solely) through the operator insertion Q{i(). The 
real part x controls the location of this insertion, 
and thereby affects both the cylinder evolution and the 
electron-quasihole commutation sign. Second, we note 
that the amount of U(l) charge carried by an Abelian 


quasihole operator Q(fC) is exactly opposite to the charge 
carried by an empty orbital = U ( 7 ). This can be 
seen by comparing the exponents in Eqs. (30) and (25). 
Therefore, to keep the CFT boundary states fixed, upon 
inserting a quasihole we need to increase the total num¬ 
ber of Landau orbitals by one. Third, in the presence 
of multiple quasiholes, the conformal correlator exhibits 
a nontrivial monodromy as a function of quasihole posi¬ 
tions. In the MPS, this monodromy manifests itself as 


k 


k + r 


(44) 


The electronic correlations in such states are character¬ 
ized by the presence of A:-particle clusters. The fermionic 
wave function is given by the product of a Jastrow fac¬ 
tor and a bosonic part. The latter bosonic wave function 
vanishes to order r when {k + 1 ) particles come together. 
The r = 2 case corresponds to the Read-Rezayi se¬ 
ries,!^ with k = 2 being the Moore-Read state,!!^ while 
the (fc,r) = (2,3) case is the so-called Gaffnian wave 
function.!^!! These wave functions can be constructed from 
chiral conformal correlators of primary fields representing 
electrons and quasiholes.!^ The electron operator takes 
the tensor product form 


V{z) =:e' (45) 

where in addition to the free-boson vertex operator, we 
also have a primary field ip in the so-called neutral CFT. 
The fundamental quasihole is represented by 

Q(,y) = a( 7 )® (46) 


with a being another primary field in the neutral CFT. 
The reduced exponent in the boson vertex operator re¬ 
flects the fact that the fundamental quasihole is a fur¬ 
ther A:-fold fractionalization of the Abelian quasihole in 

Eq. ([^. 

The MPS auxiliary space is now given by the direct 
product of the truncated Hilbert spaces of the neutral 
GET and the free boson. As noted in Sec. [Hj the free 
boson Hilbert space can be naturally broken into sectors 
labeled by the conserved U(l) charge gq. The neutral 
GET Hilbert space can be similarly split into different 
representations of the neutral Virasoro algebra.!^ Each 
representation, called a Verma module, is spanned by the 
conformal family of Virasoro descendants generated from 
a single primary state. Therefore, each Verma module 
of the neutral CFT is labeled by a primary field, which 
we refer to as the “topological charge” and denote by 
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Latin indices a,b,c,.... Taken together, the MPS auxil¬ 
iary space can be split into different sectors labeled by the 
topological and the U(l) chargesP^^For later convenience, 
we define the projector into a single Verma module c as 

P(c)=^|a)(a|. (47) 

aGc 


approach is numerically inefficient due to the complex¬ 
ity of the extended algebra and the proliferation of null 
modes. Hence, in the following we stick to the Virasoro 
description, but for succinctness we keep W and ip im¬ 
plicit whenever possible. With this caveat in mind, the 
fusion rules involving the electron 'ijj = 'ipi or the quasi¬ 
hole cr = (Ti are given by 


As shown in Sec. EH a quasihole at x + K is repre¬ 
sented by the insertion of Q{iC) into the cylinder evolu¬ 
tion at time slice y, as in Eq. (36). There are two extra 
complications for the non-Abelian case in Eq. (46). First, 
we need to resolve the topologically degenerate states 
associated with multiple pinned non-Abelian quasiholes. 
Second, we need to generalize the anti-commutativity be¬ 
tween electrons and Abelian quasiholes in Eq. (38). As 
noted earlier, a non-Abelian quasihole can be seen as a 
further fc-fold fractionalization of an Abelian quasihole. 
This hints at a fc-way split of the anti-commutation mi¬ 
nus sign. However, to have a single-valued electron wave 
function, the commutation phase between an electron 
and a quasihole must square to unity. As we demon¬ 
strate below, the solution to this conundrum turns out 
to be letting each of the k parts of a quasihole anticom¬ 
mute with only one out of every k electrons. 


A. Neutral CFT examples 

Before diving into the details of the quasihole operator, 
we first go through the field content and the fusion rules 
of the neutral CFT for a few representative theories.^^^IIIl 

Moore-Read. The neutral CFT for the Moore-Read 
Pfaffian stat^i^ is the minimal model A4(3,4) with cen¬ 
tral charge c= \. The primary fields (1,'0,cr) have scal¬ 
ing dimensions (0, j|^). The fusion rules are given by 

if X = 1, a X Ip = a, a X a = 1 Ip. (48) 

Gaffnian. The neutral CFT for the Gaffnian wave 
functiorPi^ is the nonunitary minimal model A4(3,5), 
with a negative central charge c = — |. The primary 
fields {l,ip,a,(p) have scaling dimensions ( 0 , |,—^,|). 
The fusion rules involving ip or a are given by 

iljxip = l, ipx a = ip, ipx p = a, 
a X a = 1 If, a X ip = if a. 

Z 3 Read-Rezayi. The neutral CFT for the Z 3 Read- 
Rezayi state is the Z 3 parafermionic variant!^ of the 
c = I minimal model A4(5,6). The Virasoro primary 
fields {l,ipi,ip 2 ,W,e,ai,a 2 ,p) have scaling dimensions 
(0, f, §,3, |, 3 ^, 1) and Z 3 charges (0,2,1,0,0,1,2,0). 

This theory actually enjoys an extended W 3 algebrrPSl be¬ 
yond the Virasoro. In this language, the W field is not 
a primary field, but rather a descendant of the identity 
under the larger W 3 algebra. Similarly, p is actually a 
descendant of e. However, as noted in Ref. [m the Ws 



ifl 

i >2 

e 

<Jl 

0-2 

1 pl 

V '2 

1 

0-2 

e 

CTl 

0-1 

e 

0-2 

if 2 CTl 

ifl + <72 

1 - 1 - £ 


Note that the a x a fusion has multiple outcomes 
in all of the examples above. Such fusion ambiguity 
is characteristic of the CFT representation of the non- 
Abelian quasiholes. For completeness, we list the struc¬ 
ture constants for the so-called operator product expan¬ 
sion needed to construct the MPS matrix elements in 
Appendix 


B. Conformal blocks and fusion trees 

We consider conformal correlators with multiple non- 
Abelian quasihole insertions 




(50) 


Due to the nontrivial fusion rule of the a field in Q{t]) 
[Eq. (46)], for each set of quasihole coordinates, the above 
expression does not produce a single wave function. In- 
stead, it defines a vector space of degenerate wave func¬ 
tions.^ A set of basis states in this space, called confor¬ 
mal blocks, can be obtained by specifying how the h elds 
fuse together in terms of a fusion tree diagram.l^SMl 
only need to consider the neutral CFT, since the free bo¬ 
son has trivial fusion rules dictated by the U(l) charge 
conservation. The structure of the fusion tree reflects 
the ordering of the successive fusions, and different struc¬ 
tures correspond to different basis choices for the same 
vector space. For a given fusion tree structure, the corre¬ 
sponding conformal-block basis states can be enumerated 
by finding all the topological charge labelings compatible 
with the fusion rules. 

We now consider the MPS representation of the con¬ 
formal blocks associated with Eq. (50). This mandates 


the fields to be fused sequentially in time ordering into 
the |in) state. As explained in Sec. HE each quasihole is 


pinned at a particular time slice while the electrons need 
to be placed at the center of each occupied orbital. This 
requires us to blend the quasihole operators into the chain 
of electron operators. Specifically, we have a Vo operator 
placed at the center of each occupied Landau orbital and 
a Q{iC,) operator inserted at each quasihole position y, 
and the operators are lined up in time ordering. This 
leads to a fusion tree with a linear structure. Consider 






















for example the following amplitude for the Moore-Read 
state 


•VoQ«i)VoQ«2)|1). 


(51) 


To reduce clutter, here we have omitted the interleaving 
cylinder evolution operators U{s). The fusion tree takes 
the form 


tp a 'ip a 


1 


(52) 


where the undecided topological charges (a, b) could be 
either (1,'0) or ('0,1) according to the fusion rules. It 
should be noted that the imaginary time x points in the 
left direction in the above diagram, in accordance with 
the operator time ordering. To construct the conformal 
block for either choice of the topological charges (o, b), we 
need to materialize the fusion channel choice in Eq. (511. 


This amounts to inserting a Verma module projector to 
the left of each field insertion, 

(• • • r{b) Vo r{a) Q(*Ci) Via) Vq V{a) Q« 2 ) |1). (53) 

It should be noted that the placement of the electron 
Vo operators depends on how and where the Landau or¬ 
bitals are occupied. As a result, with quasihole inser¬ 
tions in the bulk, the fusion trees for different Slater- 
determinant amplitudes naturally have different order¬ 
ings of their quasihole and electron “branches”. For ex¬ 
ample, with a non-Abelian quasihole pinned at t] = I 7 , 
the following two Slater-determinant amplitudes have dif¬ 
ferent operator ordering and thereby different fusion tree 
structure, 


■ 101010,11 
■01101101 


•VoVoVoQ(O)VoVoll), 

•VoVoVoVoQ(O)Voll). 


(54) 


C. Exchanging two branches 


To connect the fusion tree with the desired structure 


in Eq. (55) to the primitive time-ordered tree in Eq. (521, 


we need to move the quasihole insertions across electron 
operators into the bulk. The elementary move is 


iP 


a 


if 


d 


d 


(56) 



These two fusion trees span the same vector space, and 
they can be related by a linear transform 


(57) 


which is nothing but the CFT half-braid matrix between 
the electron and the quasihole fields. It should be noted 
that the U(I) part of the CFT also contributes a phase 
factor, even though it has trivial fusions and has been 
omitted from the above diagrams. Also, for the CFT cor¬ 
relator to represent a physical wave function, the electron 
operator must be local with respect to the quasiholes. 
As a result, it does not matter whether the electron- 
quasihole half braid is done clockwise or counterclock¬ 
wise. This brings Eq. (57) to a more familiar form. 


-0 

= E [b: 


■ <7'lpC' 



(58) 


Formally, the half-braid matrix can be decom¬ 

posed into the so-called F and R moves of the direct- 
product CFT.I^ Again, the contribution from the U(I) 
part must be included despite its omission from the di¬ 
agrams. The fusion F matrix is a generalization of the 
Wigner 6j-symbol, 


Here the occupation numbers are listed in reverse to be 
consistent with the time ordering, the star marks the po¬ 
sition of the quasihole relative to the Landau orbitals, 
and again we have omitted the interleaving cylinder evo¬ 
lutions in order to highlight the operator ordering. Now, 
to obtain the correct wave function, we must resolve 
every Slater-determinant amplitude in the same fusion 
tree basis. A natural choice is to have all the quasihole 
fields placed at the beginning (rightmost end) of the linear 
strueture. For the example given above, we want 



For each Slater-determinant amplitude, we need to reor¬ 
ganize the time-ordered fusion tree in Eq. (52), by bring¬ 
ing all the quasihole branches across the electrons to the 
rightmost, while keeping their relative ordering. As we 
show next, this reshuffling transformation amounts to 
adding a particular minus sign to each quasihole oper¬ 
ator, generalizing the prescription for the Abelian quasi¬ 
hole in Eq. (43). 


(59) 

li ® » d 

where g is summed over all topological charges compat¬ 
ible with the fusion rules. And the R matrix gives the 
exchange phase in a definite fusion channel, 

a b a b 

y = Y ■ 

9 9 

Note that the fusion tree layout in the above definitions 
is slightly different from the standard convention in the 
literature.l^ This is due to our choice of pointing the 
imaginary time x in the left direction in accordance with 
the operator time ordering. Composing the F and R 
moves, we find 



(61) 
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The transformation coefficient in the above equation is 
nothing but In principle, we can solve the pen¬ 

tagon and the hexagon equation^^ for the F and the R 
matrices, and fix the gauge according to the structure 
constants in the operator product expansion of the chiral 
fields (see Appendix]^. This would be a time-consuming 
task. In practice, however, we can more easily determine 
the matrix from a simple numerical calculation. 

To this end, we consider the conformal-block version of 
Eq. (|57l). 


{d\V{iy)V{e)Q{iO\c) = 

E [BT]e}(d\Q{^OV{^)V{^y)\c). (62) 

/ 


Here, we place the operators at adjacent time slices to 
avoid complications from the cylinder evolution, and for 
simplicity we choose the end states (d] and |c) to be the 
primary states in the corresponding Verma modules. Re¬ 
call that the core part of the MPS implementation is 
a truncated CFT calculator. By numerically evaluating 
the two (truncated) correlators in the above equation as 
functions of {y, (C), we can extract the matrix in the 
gauge determined by the chiral structure constants. 

We first consider the Moore-Read and the Gaffnian 
states. They are clustered wave functions with k = 2. 
For these theories, the fusion tpxa always yields a unique 
result for any topological charge a. Therefore, the topo¬ 
logical charges d and / in Eq. (571 are completely fixed 
by the choice of c and e, namely 


d = ijj X e, f = Ip ^ c, (63) 

with e among the (possibly) multiple fusion outcomes of 
axe. This allows us to adopt the shorthand notation 

[BT]ef ^ Bl, (64) 

and it reduces Eq. ( [^ to 

(V^ X e| V(iy) V{e) Q«) |c) = 

Bl ■ {pj X e| Q{iC) Viip X c) V{iy) |c). (65) 


The single-valuedness of the electron wave function dic¬ 
tates that Bg has to be ±1. For the Moore-Read state, 
from the numerics we find 


as primary fields, even the electron field ip = ipi has 
nontrivial fusion rules. 


tpi X 'tp 2 = 1 + W, -ipi X ai = e + (f, ( 68 ) 

which would not bode well for our proposed method. 
However, the implicit W 3 symmetry forbids us to split 
the channels on the right-hand side into separate con¬ 
formal blocks. To take advantage of this, we bind the 
Verma module IT to 1 and bind to e by redefining the 
projectors [Eq. HI)], 

E E (69) 

aGl,W OL^s,!^ 


and then we simply forget about W and (p when labeling 
the fusion trees. Now that W and ip are formally gone, 
we are allowed to use the shorthand notation in Eq. (641. 


Crucially, we verify from the numerics that after such 
patching (but not before), the reduced linear transform 
in Eq. (65) still holds. This compatibility can be at¬ 


tributed to the underlying extended W 3 algebra. The B^ 
coefficients are found to be 


Bl^ = Bp = Bp = Bp = Bp = Bp = 1, 
Bp = Bp = Bp = - 1 . 


(70) 


D. Reshuffling quasiholes into the bulk 


Through successive applications of the elementary 
move in Eq. (57), we can achieve the reshuffling 



across an arbitrary number n of electrons. Note that the 
topological charges at the top and the bottom must agree 
between the two trees; otherwise the conformal blocks be¬ 
long to distinct linear spaces. Also, due to the trivial fu¬ 
sion rules of the ip field in the Moore-Read, the Gaffnian, 
and the Z 3 Read-Rezayi states, all the topological charges 
are completely fixed by c and e. In particular, we deduce 
that the in-situ fusion channels for the cr field after the 
reshuffling are given by the successive fusion with n ip 
fields. 


e' = ip'^ X e, c’ = Ip"' X c, (72) 


Bi=B5 = l, Bt=Bl = -l. (66) 

And for the Gaffnian state, we have 

B1=B; = B- = B^ = 1, = = (67) 

Depending on the fusion channel context, the electron 
and the quasihole fields either commute or anticommute. 

For the Zfc ^3 Read-Rezayi state, due to our choice of 
keeping the W 3 algebra implicit and treating W and ip 


and the corresponding projected quasihole operator is 

Q{e,c,np() = V(ip" x e) Q«) V{ip" x c). (73) 

The linear transform between the two trees in Eq. 0 
consists of a simple sign factor. 












10 


Here, only a subset of the n V' fields anticommute with 
a, and 


N{e, c,n) G Z 


(75) 


counts the size of this anticommuting subset. Consider 
(e, c, n) = (cr, 1,4) for the Moore-Read state as an exam¬ 
ple. Recall from Eq. ( 66 ) that = 1 and = —1. 
Using Eq. (57), we find 


Ip ip ip Ip 



Ip Ip Ip 

^ = Bl\ - 


ip Ip 

= BlBt\ 




= BPB^BPB. 


In the above commuting process, the a held picks up 
a minus sign Bf = — 1 from only every other ip field 
(dashed line). This is a direct consequence of Eq. ( 66 ). 


The alternating pattern in the chain of ip fields has peri¬ 
odicity two, as expected from the fusion rule ip'^ = 1. 

Eollowing the above procedure, for each channel choice 
(e, c), we can identify the subset of ip fields that anticom¬ 
mute with a (dashed lines). For the Moore-Read state 
we find. 


Ip Ip Ip Ip a 



Ip Ip Ip Ip 



(77) 


Counting the dashed lines, we have for Eq. (74) 
N{a, 1, n) = N{ip, cr, n) = [§J , 

7V(cr, Ip, n) = N{l,a, n) = , 


(78) 



Ip Ip Ip Ip a 


where the floor function [tj G Z denotes the largest in¬ 
teger no greater than t G R. Similarly, for the Gaffnian 
state we have 

Ip Ip Ip Ip a 
a 


(79) 


Ip Ip Ip Ip cr Ip Ip Ip Ip a 

And the number of anticommuting ip fields is given by 
N{a, l,n) = N{ip, tp, n) = [f J , 

N{ip,ip,n) = N{l,a,n) = , 

N{ip,a,n) =N{a,(p,n) =0. 



Note that for (e,c) = (</ 5 ,cr) and {<T,(p), the cr field com¬ 
mute with all Ip fields without any minus sign. For the 
Z 3 Read-Rezayi state we find. 


1 pl 1pl 1pl 1pl 1pl 1pi_ (Ti 



C 


V>1 Ipl Ipl Ipl Ipl Ipl (Ti 



C 


Ipl Ipl Ipl Ipl Ipl Ipl (7l 



C 


c 


1 


CTl 



e 


.0-1 


.i’l. 


e 

c 


Ipl 


e 


CTl 

e 


e 


.■02. 


.^2. 

c 


V'2 


^2 


e 

e 


.CT2. 


1 


.CTl 


( 81 ) 

And the number of anti-commuting ipi fields is given by 


N(ai, l,n) = N(ipi,ai,n) = N{e, era, n) = [|J , 

N{e, ipi,n) = N{ip 2 , e, n) = N{a 2 , cri, n) = , (82) 

N{a 2 , ip 2 , n) = iV(l, 0 - 2 , n) = fV(cri, e, n) = . 

Again, we emphasize that the sign structure has contri¬ 
butions from both the neutral and the implicit U(l) parts 
of the direct-product CFT. 


E. Putting the pieces together 


We apply the reshuffling formula (74) to bring the 


quasihole and the electron fields into time ordering. A 
quasihole at position x + *C should be placed between the 
orbitals = [x/xj and + 1 [Eq. (32)]. Starting from 
the rightmost end of the fusion tree, to reach this time- 
ordered position, the number of electrons it needs to cross 
is given by the number of occupied orbitals with center 
X < X- Consider m fundamental quasiholes with order¬ 
ing Xi > ■ ■ ■ > Xm [Eq. (41)]. The number of occupied 
orbitals with center a; < Xi is counted by the operator 


hi = y/vao + Uxi + 1 )J^ - 


m — I 
k ' 


(83) 


inserted between the C’" operators for orbitals and 
jy^ + 1. The above formula is adapted from the Abelian 
case [Eq. mi)]- The extra k in the denominator of the 
last term reflects the further /c-fold fractionalization of 
the U(l) charge of a fundamental quasihole in the non- 
Abelian states [Eq. ([4^]. 

Now we are finally ready to synthesize the full expres¬ 
sion for the non-Abelian quasihole operator. We would 
like to obtain the MPS description of the conformal block 
with m fundamental quasiholes specified by the fusion 
tree 

rn 



(80) 


( 84 ) 
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(a) 



(b) 

e + iC - -e + i( 

f- *- 


FIG. 2. (a) Braiding one quasihole around another. The 

bystander quasiholes in the shaded regions do not participate 
in the braiding process other than setting the fusion channel 
context. The vertical dotted line marks the location of the 
branch cut discontinuity, (b) An infinitesimal segment of the 
braiding path crossing the branch cut from above. 


Here the quasihole fields are placed at the rightmost end 
of the linear tree structure as in Eq. (55), and we have 


omitted the electron branches on the left side. The topo¬ 
logical charges satisfy bi-i G a x bi, with 6 ^ = 1. In the 
MPS, the quasiholes are placed in time ordering, and the 
Z-th quasihole with fusion context ( 6 ;_i, bi) is represented 
by the insertion of 

U{5xi) Q{bi-i,bi,hi, iCi) mSxi)]-^ A) 

(85) 

between the C"" operators for orbitals and -h 1. In 
the above expression, the displacement Sxi is defined in 
Eq. (34), the projected quasihole operator Q(e,c,n,iC) 


is defined in Eq. (73), the occupied orbital counter fii 


is defined in Eq. (83), and the anticommuting subset 


counter N{e,c,n) defined in Eq. (74) has values given 
by Eqs. ([78|), (|l0|, and (|82|). The above MPS representa¬ 


tion preserves the conformal-block normalization for the 
non-Abelian quasiholes, and it is the main result of this 
paper. 


F. Branch cut structure 

As the final note, in this section we examine the mon- 
odromy structure of the MPS when we braid one non- 
Abelian quasihole around another (Fig. [^. The MPS 
prescription produces conformal blocks in a particular 
fusion tree basis, where the quasiholes are placed in time 
ordering at the rightmost end of the linear structure be¬ 
fore the electrons. This allows us to focus on the fusion 
tree segment actively involved in the braiding process 


a a 



without worrying about the bystander quasiholes (Fig.[^ 
or the electrons, which just define the topological charges 
(a, c). Due to the time-ordered insertion of the quasi¬ 
holes, each conformal block develops a branch cut discon¬ 
tinuity when the two quasiholes coincide in the horizontal 
direction, as illustrated by the dotted lines in Fig. [^a). 
This discontinuity represents an abrupt change of basis 


and does not correspond to a physical singularity. To 
understand how the conformal blocks across the branch 
cut are related, we consider the two points infinitesimally 
close to the cut, as shown in Fig. [^b). To be specific, 
we place the stationary quasihole at ryo = 0 and examine 
the infinitesimal segment across the cut above it, from 
—e + if to e -I- if. The conformal blocks on the two sides 
of the cut are given by 


a{e+if) (t(0) (j(0) c7{-e+if) 



In the limit of e —)■ 0+, the conformal blocks on the right 
side can be related via continuity to the twisted trees 


cr( 0 ) a{—e+if) 

"a 


(T(e-NC) o'(O) 



( 88 ) 


Further, we can untwist the fusion trees using the half¬ 
braid B matrix of the direct-product CFT [Eq. (58)], 


a{e+if) cr( 0 ) 

A = E [«” 



a{e+if) a( 0 ) 


J bd 


(89) 


leading to 

(t(0) a{-e+iC) 



a(e+if) a( 0 ) 



Similarly, across the branch cut below the stationary 
quasihole we find 


cr(e-iC) o-(O) 



E 


cr( 0 ) a(-e-iC) 



Therefore, the MPS conformal blocks across the branch 
cuts are related by the half-braid B matrices of the direct- 
product CFT. In other words, the full-braid monodromy 
of the conformal blocks are concentrated in the branch 
cut singularities. This is a special feature of our fusion 
tree basis with time ordering, and it drastically simplifies 
the microscopic demonstration of the topological nature 
of the quasihole braiding statistics.^ With the branch 
cuts contributing the CFT monodromy matrix, now we 
just need to show that the braiding process away from the 
branch cuts accumulates only a simple Aharonov-Bohm 
phase. For the quasihole wave functions constructed from 
the conformal correlators, the latter condition can be fur¬ 
ther reduced to the exponential convergence of the over¬ 
lap matrix between different conformal blocks at large 
quasihole separations.^^ In an earlier paper, we took ad¬ 
vantage of this line of simplifications to demonstrate the 
Fibonacci nature of the Z 3 Read-Rezayi quasiholes. This 
exploitation depends crucially on the fact that our MPS 
construction is a literal transcription of the conformal 
blocks preserving both the monodromy structure and the 
plasma normalization. 
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IV. CONCLUSION 


into Eq. (A3), we find 


In this paper we have presented a pedagogical deriva¬ 
tion of the conformal-block MPS for non-Abelian quasi¬ 
holes. The procedure is exemplified using the Moore- 
Read, the Gaffnian, and the Z3 Read-Rezayi states. Our 
prescription preserves the monodromy structure and the 
plasma normalization of the conformal blocks, and the re¬ 
sulting MPS explicitly manifests the putative quasihole 
braiding statistics as half-monodromy matrices across 
branch cuts. 
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Appendix A: Uniform background charge 


In this Appendix we prove Eq. (12), namely that the 
inclusion of the uniform background charge [Eq. 0] 


C>bc = : exp 



(AI) 


in the conformal correlator correctly produces the 
Landau-gauge Gaussian factor with an extra non- 
holomorphic gauge transform, such that 


fix,y) = 


t 


dx'dy' log 


2 . ^j{x + iy — x' — iy') 
7 2 


(A5) 


The integration in the above equation is performed over 
the cylinder surface {x',y') G K x {—^Ly, \Ly). We pro¬ 
ceed by separating the real and the imaginary parts of 
the complex logarithm, logit = log |it| -I- ilmlogit. 

The real part of the logarithm is nothing but the 
Coulomb potential on the cylinder geometry,123 in the 
sense that 


VMog 


— sinh 
7 


l{z - z') 
2 


27r S(z — z'), 


(A6) 


with cylinder identification z z + iLy. The real part 
of f{x,y) is thus the scalar potential due to the uniform 
background charge, satisfying 


V^Ref{x,y) = - J dx'dy' 6{x-x')6{y-y') = -1. (A7) 


Imposing the reflection and the rotational symmetries on 
the cylinder, 

^ef{x, y) = Re/(-a:, y), dyRef{x, y) = 0, (A8) 

we can integrate the above Poisson equation and obtain 

Re/(a;,y) =-y. (A9) 


^(zi,..., z„) = (V(zi) • • • V(z„) Obc) (A2) 

is a legitimate many-body wave function in the lowest 
Landau level with the Landau gauge. 

The above statement applies to quantum Hall states 
described by a generic direct-product GET. Without loss 
of generality, we can focus on the U(I) part of the electron 

operator, V{z) :, since the background charge 

Obc does not couple to the neutral GET. Further, due to 
the noninteracting nature of the U(l) boson, we only need 
to consider the contraction of a single electron operator 
with the background charge. 


As claimed, this correctly reproduces the Gaussian factor 
in the one-body Landau orbital [Eq. (10)]. 

For the imaginary part of the logarithm, we employ 
the identitji^ 


T . . ^ 1 [x + iy) y + nLy 

Im log sinh - -^^ arctan - - -. (AIO) 

n 

It is not hard to see why this holds: starting from 
Im log '^{x + iy) = arctan ^, the periodic structure intro¬ 
duced by the hyperbolic sine is matched by the periodic 
sum over n. As a result, the imaginary part of f{x,y) is 
given by 


g/G.y) 




= I da;'dy' rp{x'+iy') . 


dx'dy' {(j){x + iy)4>ix' + iy')) j ■ 

(A3) 

We work on the cylinder geometry with perimeter Ly 
and inverse radius 7 = 27:/Ly. Plugging the boson prop¬ 
agator 





log 


2 sinh ^ 
7 2 


(A4) 


Im/(a;,y) = —y f dx'dy'^ arctan^ + nLy 


cyl 


— / dx' [ dy' arctan -—— 

27r Jm Jr X- x' 

— ^ / dx'sign(x — x') = —xy. 

^ .m 


X — X' 
7 


Here, we have joined the cylinder integrals in the infinite 
sum to tile the K x K plane, with a symmetric regular¬ 
ization for the ±00 limits. 
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To sum up, we have shown that the contraction with 
the cylinder background charge is given by 


• e* dx'dy' cl>{x' +iv') . 


= (A12) 


This proves that the wave function defined in Eq. (12) in¬ 


deed lives in the Landau gauge with the correct Gaussian 
factor. 


Appendix B: Calculating matrix elements 


As noted in Sec. IID on a cylinder with finite perime¬ 
ter Ly = 27r/7 the evolution factor e”'*' allows us to 
truncate the conformal Hilbert space according to Lq 
eigenvalues. In the actual calculations behind our previ¬ 
ous paper,we kept the descendant states with Lq < 14 
for the Moore-Read and the Gaffnian states, and Lq < 13 
for the Z3 Read-Rezayi state. This is necessary to reach 
convergence for cylinder perimeter Ly up to 25 magnetic 
lengths. The resulting MPS auxiliary space has size up 
to 3 X 10^. The calculation of physical observables is 
carried out over the direct product of two copies of the 
auxiliary space, one for (('I'| and one for |'I')). The size of 
this direct-product space can be close to 10®. 


In the main text the MPS is described in terms of oper¬ 
ators acting on the GET Hilbert space. In this Appendix, 
we briefly describe the construction of the matrix repre¬ 
sentation of these operators. A detailed procedure for 
the calculation of the primary field matrix elements has 
been published in Ref. [TTJ Leveraging this result, all we 
need here is to map the operators of interest from the 
cylinder geometry to the plane. 

Under the conformal map 2: —> from the cylinder 

to the plane, a primary field 4> with scaling dimension 
A$ transforms by 

$(2) = $pia„e(e^"). (HI) 


Appendix C: Operator product expansions 

In the main text we have listed the field content and 
the fusion rules for the Moore-Read state, the Z3 Read- 
Rezayi state, and the Gaffnian wave function. To actu¬ 
ally construct the MPS matrices, we also need the full 
operator product expansion with structure constants as 
noted in the previous appendix. In the following, we use 
a shorthand notation 

X ^ (Cl) 

I 


Between two generic GET energy eigenstates |a) and |/3) 
with Lq eigenvalues Aq and A^g, the matrix element of 
^{iy) reads 


for the operator product expansion 

I 


{amtym = 

= e*^^(^“-A^)7^-(a|<i>pia„e(l)|/3). 

The planar matrix element (ct|4>piane(l)l/3) can be further 
related through algebraic manipulationsIH to the struc¬ 
ture constant associated with the $q, terms in the opera¬ 
tor product expansion 4) x where <I>q, and are the 
parent primary fields for the states |a) and \/3). This pro¬ 
cedure can be directly applied to the quasihole Q(iC)- As 
for the electron operator, we can reduce the zero mode 
Vo [Eq. @] to a similar form, 

(a|Vo|^) =7^" ^ 

Jo 

= 7^<5A„,Ap(a|Vplane(l)|/3)- 

(B3) 

The relevant structure constants of the neutral GET for 
various quantum Hall states are documented in the next 
appendix. 


where hi is the scaling dimension of the primary field 
Eor the Moore-Read state, we have 


tp X tp = 1, a X Ip = —j= <J, a X a = 1 H— -j= xp. 

v2 V 2 

Eor the Gaffnian state, we have 

ipxpj = l, tpxa = ^ip, ipx(p=^a, 

a X a = 1 + Cl (fi, a X if = —^ip + Ci a, 

V2 

with the constant Ci given by 

r = / %/5-l ^( 5 ) 

' ' ' ^r(i)r(i) 


(C3) 


(C4) 


(C5) 


Einally, for the Z3 Read-Rezayi state, the structure con¬ 
stants are listed in Table HI 


^ D. C. Tsui, H. L. Stormer, and A. C. Gossard, Phys. Rev. ^ X.-G. Wen, Advances in Physics 44, 405 (1995) 
Lett. 48, 1559 (1982) 
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TABLE I. The OPE structure constants for the Z 3 Read- 
Rezayi GET. Listed in column a, row b is th e fusion result 
axb using the shorthand notation in Eq. (Cl I. The constant 
C 2 is given 


im 

rr(f)l 

i/r(|) 

k(§)J 



tpi 


Ipl 


\/i ^ + l\/i F 

■02 

1-^W 


w 



e 

\/i '^2 

\J\ V’2 + \/^ O' 1 

CJ-1 


^'01 + V 2 G 2 02 

CJ-2 



T 

1 ^ 

3 Y 2 ^2 

“l\/i ^>2 + 
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